Conditions for Vanishing Central-well Population in Triple-well Adiabatic Transport 



OO 
O 

o 

(N 

> 

o 
O 

ON 



c 

o 
o 



> 

(N 

m 
m 



oo 
o 



X 



Tomas Opatrny 1 and Kunal K. Das 2 ' 3 
1 Optics Department, Faculty of Science, Palacky University, 
17. Listopadu 50, 77200 Olomouc, Czech Republic 
2 Department of Physical Sciences, Kutztown University of Pennsylvania, Kutztown, PA 19530, USA and 
3 Department of Physics, Fordham University, Bronx, New York 10458, USA 
(Dated: October 19, 2008) 

Analytical expressions are derived for coherent tunneling via adiabatic passage (CTAP) in a triple 
well system with negligible central-well population at all times during the transfer. It is shown that 
a manipulation of the depths of the extreme- wells, correlated with the time variation of the non- 
adjacent barriers is essential for maintaining vanishing population of the central well. The validity 
of our conditions are demonstrated with a numerical solution of the time-dependent Schrodinger 
equation. The transfer process is interpreted in terms of a current through the central well. 

PACS numbers: 05.60.Gg, 03.65.Xp,73.63.Hs, 73.23.-b 



I. INTRODUCTION 

Certain characters in the Harry Potter adventures have 
the magical power to move from one room to another 
without having to go through the staircase in between 
fl[ . Even without resorting to magic and wizardry, some- 
thing quite similar can be brought to pass in the quan- 
tum world: among three adjacent potential wells a, b, 
and c (see Fig. [T]), a particle in well a at one end can be 
transferred to the non- adjacent well c at the other end 
without bypassing or occupying the well b, that lies in 
between. The principle behind such a transfer has been 
well established for quite some time in the context of 
internal atomic states for a two-photon process known 
as STIRAP (stimulated Raman adiabatic pumping) [2], 
wherein atoms are pumped from a state a in their elec- 
tronic ground-state manifold to another state c in the 
same manifold by means of a pair of time varying laser 
fields that couple a and c to an excited state b. 

The idea can be applied to other similar manifesta- 
tions of discrete quantum states; specifically when those 
states represent the spatial wavefunctions of physical par- 
ticles, one can achieve the seemingly magical transport 
between two locations without traversing the interven- 
ing space. The last few years have witnessed a surge 
of theoretical interest in applying the STIRAP princi- 
ple to spatial transport of physical particles in a variety 
of systems: atoms in optical microtraps electrons in 
quantum dots 0, [5| , Bose- Einstein condensates [f| 0] , su- 
perconductors |8|, spin transport in the context of quan- 
tum information [9(. Experimental demonstration was 
achieved for photons in optical waveguides [To| . Even a 
new acronym CTAP (coherent tunneling via Adiabatic 
Passage) has been coined to describe the spatial ana- 
log The popular theoretical description of the pro- 
cess is borrowed from quantum optics, and is in terms 
of a discrete 3x3 Hamiltonian matrix with variable off- 
diagonal elements corresponding to the coupling between 
the traps. Here, we adopt a more flexible approach, ap- 
propriate for spatial transfer, where the evolution is de- 
scribed by a time dependent Schrodinger equation for 



the appropriate potentials forming the wells and the bar- 
riers. The advantage of such a method lies in its capa- 
bility to model arbitrary potential shapes and variations 
accurately d, 0]- 

Regardless of the method used, results on such adia- 
batic transfer have been numerical in nature, with the 
tradeoff that the parametric dependencies of the trans- 
fer process are hard to identify. Our goal in this paper 
is to obtain simple analytical conditions for the system 
parameters by considering a somewhat simplified triple- 
well model which captures the essential features of more 
complex potentials. We find that in order to have a 
vanishing population of the central well at all times cer- 
tain supplementary conditions need to be satisfied: First, 
out of the three lowest-energy and relevant eigenstates of 
the Hamiltonian, the middle-energy one should maintain 
as low a population of well b as possible, which can be 
achieved by requiring its node to be at the center of well 
b. Secondly the energy eigenvalues of the instantaneous 
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FIG. 1: Potential structure for the three- well system with 
schematic plots of the three lowest energy eigenstates. The 
three wells a, b, and c of length L are separated by two delta- 
function barriers of strengths Wi,2- The potential of the mid- 
dle well is set to 0, the outer wells a and c have variable 
potentials U a , c - 
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Hamiltonian need to be manipulated in such a way that 
the process stays in the adiabatic regime, which can be 
achieved by keeping the energy of the middle state as 
far from the adjacent states as possible. We show that 
the implementation of these conditions necessitate the 
manipulation of the well depths as well as the barrier 
heights. 

The paper is structured as follows: In Sec. II, we de- 
rive the conditions for keeping the central well popula- 
tion negligible at all times in the limiting case of delta- 
function potential barriers, and in Sec. Ill we show that 
those conditions apply for finite width barriers that are 
narrow and high at all times, and also qualitatively for 
wider barriers. These conditions are applied in a numer- 
ical simulation based on a time-dependent Schrodinger 
equation in Sec. IV to demonstrate the validity of our 
conditions. In Sec. VI, we provide an interpretation of 
the transfer process in terms of a high velocity current, 
and summarize our results in the conclusions section. 



II. SYSTEM PARAMETERS AND 
EIGENSTATES 

In STIRAP the transfer is achieved by using a coun- 
terintuitive coupling sequence, where b is first coupled to 
the empty target state c, and only when that coupling is 
reduced, the ab coupling is ramped up. Appropriate adi- 
abatic variation of the two couplings map the state b to a 
dark state superposition of states a and c, which contains 
no contribution of state b. The name 'dark state' refers to 
the fact that atoms in this superposition cannot interact 
with the laser fields. By gradually switching off the ab 
coupling, the dark state is adiabatically transformed into 
c. During the process, there is virtually no occupation 
of the intermediate state b which would spontaneously 
decay into one of lower energy states. 

In the spatial analog, the three states are imple- 
mented by three adjacent potential wells with the inverse 
strengths of the two intra- well barriers signifying the cou- 
pling between the wells. We implement this process by 
considering a one dimensional system with the potential 
energy structure as follows (see Fig. [TJ: 
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where each well spans a length L in the intervals 
a = ,b= [-\L;\L] ,c = [§£; §L] . Effec- 

tive ID behavior can be implemented for atoms in micro- 
magnetic traps fllj ] with tight transverse confinement, 



and for electrons in size-quantized nanowires 
delta-potentials are idealizations of narrow and high po- 
tential barriers, which can be implemented with tightly 
focused lasers for ultracold atoms, or sharp field-gradient 
between adjacent quantum dots in the case of electrons. 



The potential energy shifts U n r are analogs of detuning 
among internal atomic states [131 ]. Denoting the funda- 
mental wavevector for an infinite square well of length L 
by fco = tt/L for a state with wavelength A = 2L, we de- 
fine the dimcnsionlcss coupling parameters ai,2 for the 
inverse-strengths of delta-function potentials, and u a , c 
for the constant energy shifts of the two outer wells: 
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Consistent with the assumption of high barriers, we as- 
sume a 12 1 during the whole process so that only the 
lowest orders of a 1,2 need to be considered. 

The initial conditions, along with the adiabaticity, can 
ensure that only the three lowest energy states of the in- 
stantaneous Hamiltonian are energetically accessible and 
relevant; we denote the states by ip g (ground), ip d (dark), 
and ip e (excited) , in order of increasing energy. We study 
their properties below to extract conditions on the sys- 
tem parameters. 



A. Dark state 

The middle-energy eigenstate plays the crucial role in 
the CTAP time evolution, by evolving as a dark state 
with little contribution from the central well, b. In or- 
der to minimize the population of well 6, we assume a 
symmetric form of the middle-energy eigenstate 
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^Pd(x) = { B^/f sin [fe (l - y b )x] , x e b (3) 

cos [ko(l - y c )x + ^-y c ] , xec 
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The small parameters y a ,b,c <C 1 correspond to phase 
shifts of the infinite well wavevector due to the finite 
heights and widths of the intermediate barriers. The con- 
ditions of continuity of ip d at x = ±L/2, on retaining only 
terms linear in y a ,b,c lead to the relations 



A d ny a = B d = -C d iry c . 



(4) 



The time-independent Schrodinger equation applied at 
the positions of the delta walls leads to the requirements 



Va,c = "1,2, 

and applied within the wells, gives 
2mE 



(5) 
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(1 - y a y + u a = (1 - y b ) 2 = (1 - y c y + u c , (6) 



from which it follows that 

Va - tt = yb = y c 
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up to first order in y a , b ,c- This leads to the condition con- 
necting the potentials of the walls to the inverse barrier 
strengths: 



u c — u a = 2(ai - 02). 



(8) 



This relation underscores the fact that slight shifts in the 
potential wells are essential to minimize the population of 
the central well, and the optimal shifts need to be varied 
synchronously with the barrier potentials between the 
wells. If this condition is not satisfied then the dark state 
does not have the node in the middle of the central well, 
instead the node is closer to (or inside of) one or the other 
of the two walls, which due the heightened asymmetry, 
increases the net probability of the particle to be in the 
central well. Such, for instance would be the case in the 
model used in Ref. Q where the extreme wells have equal 
depths, which in our notation would mean u c = u a . 

The normalization coefficients in ([3]) are determined by 
Eqs. (0} and © to be 
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on keeping the lowest order terms in a\. 2 , and the nor- 



malization I Ad 1 2 + \Bd 



\C d \ 



1 is satisfied up to 



the same order. This assumes that the barrier heights 
are always kept high, forcing the parameters ai,2 to be 
always small. In that case, the population of the mid- 
dle well, \Bd\ 2 can be made arbitrary small compared to 
the other two wells, since the numerator of B is bilinear 
in the small parameters, but the numerators of Ad and 
Cd are linear, while they all share the same denominator 
which is effectively of linear order in the small parame- 
ters. It is clear from the expressions that by changing 
the ratio of the magnitudes of the barriers and therefore 
U\ : U2 = 012 '■ ol\ the population ratio of wells a and c 
can vary as the squares, a 2 . : a 2 , from 1 : to : 1. 



B. Ground and excited states 

The condition Eq. ([8|) fixes an optimal difference for 
the extreme well depths, but leaves the actual values free. 
Further optimization will require a knowledge of the re- 
maining states. The ground (g) and the second excited 
(e) states of the Hamiltonian both have the general form 



ip g ,e{x) = < BJ f cos [fc (l - y b )x + > 



A y I cos [ fc o(! - Va)x - ^-y a ] , x G a 

x G 6(10) 

CJ£ cos [&o(l _ Vc)x + ^-y c ] , x e c 

where we assume y a ,b,c <C 1 and (f> -C 1. The continuity 
of ipg ye at x = ±L/2 leads to the conditions 

Any a = -B (^y b + ^ ; Cny a = -B (|y b - .(11) 



The stationary Schrodinger equation applied at the po- 
sition of the delta walls leads to the requirements 
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while energy eigenvalues in the wells lead to a condition 
identical to Eq. ([7]) for the dark state. Equations (fTTj) 
yield the following relations 



Ay a), 

Ay a + By b + Cy c = 0. 



(13a) 
(13b) 



Thus all the parameters of the wavefunction are deter- 
mined by the shift y b of the wavenumber in the middle 
well, which we now solve for after simplifying the nota- 
tion by dropping the subscript: y = y b . Expressing y a c 
by means of Eq. 0, and A, C by means of Eqs. fTT j) -([T2] ) 
we find the requirement for y to be 
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Imposing Eq. ([8]) for minimum middle well population, 
the potential shifts of the well can be written as 

u a = u + 2ai; u c ^u + 2a 2 . (15) 

allowing Eq. fT4")) to be rewritten as a cubic equation 
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(qi + a.2) + a\ + a\ 
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provided that y 7^ — ^ , where an equality would lead to 
zero denominators in Eq. (fT4")) . As can be easily checked, 
y = —Mo/2 is also a root of Eq. (fT6|) , which corresponds 
to the value of y b of the dark state tpd- Therefore the 
cubic equation can be reduced to a quadratic equation 
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+ a 2 + a|,(16) 



where the value y + corresponds to ip g and y_ to ip e 



C. Choice of the well depths 

The ground and the excited states satisfying the opti- 
mal condition Eq. ((SJ) have now been determined in terms 
of the parameter uq which which remains a free parame- 
ter, and therefore allows for further optimization. That 
will be based on the requirement that the process stays as 
close to adiabatic as possible, with negligible transition 
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from the dark state to the other two available states. A 
measure of the transition rate is given by the parameter 



|<v<|ff|V<> - (MHHdW 



(17) 



where energy gap between dark state (d) and the adjacent 
states (i = g, e) sets the scale @. 

A rigorous optimization of uq with respect to (|T7|) can 
be done numerically, but our goal here is to get intuitive 
analytical expressions for the system parameters. Since a 
smaller value of A(ipi,il>d) is better in the sense of lesser 
transfer to the other states, we can deduce simple re- 
lations by requiring the denominators in (|17p for both 
tp g and ip e to be as large as possible. For A(ipi, ipd) to 
be small for both the states, the energy of the dark state 
should be exactly between the energies of the ground and 
excited states: 

4[(i-z/+) 2 + (i-y-) 21 
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(18) 



from which it follows that (up to the first order in a.\ 2) 
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~2(ai + a 2 ), 



(19) 



so that the optimal choices of potentials u a and u c are 
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(20) 



These relations show that to maintain minimum popula- 
tion in the central well, one needs to vary the well depths 
of the extreme wells in inverse proportion to the barrier 
heights. But, remarkably, the variation of the depth of 
each extreme well has to be synchronized with the height 
of the inter-well barrier which is non-adjacent to it. 
These conditions also fix y± 



y± = a\ + a 2 ± 



l 2- 



(21) 



which in turn determine the instantaneous energies of the 
three states: 
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(22b) 



as well as the parameters that determine the instanta- 
neous forms of the ground and the excited states, respec- 
tively, 
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The + corresponds to the ground state and — to the 
excited state. The knowledge of these states and their 
energies are essential for optimal time evolution. 



III. FINITE WIDTH BARRIERS 

We now show that the delta function potentials capture 
the behavior of high and narrow finite width inter-well 
barriers quite well, and the conclusions derived above 
remain qualitatively unchanged even with wider barriers. 
We relax the delta potentials to consider square barriers 
of width D, and height Ui }2 and as before introduce small 
dimensionless parameters a 1,2, now defined to be 
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We search for the conditions that would allow for the dark 
state of the sine form as in the previous case, requiring 
the wavefunction to have the form 
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D\. The same coefficient B is chosen for the barriers 
and the middle well, because up to the second order the 
derivatives are zero at the boundaries. The equivalence 
of the energy of the stationary states in all the wells leads 
to a condition similar to Eq. ([5]) for the delta barriers 



(25) 



The only difference is that in order to allow for finite 
width, we retain second order terms in y a}C , assuming 
that they are still relatively small. The continuity condi- 
tions for the wavefunction and its derivative leads to the 
equations for y a c : 



Ki }2 tanh (ki, 2 D) 

Kl,2 



kg (I - y a , c ) 
tan[7r(l - y a . c )} ' 



with — — — w [jrkoDai^] 1 — 1. 



(26) 



For narrow (k D <C 1) and high (Z7 1>2 ^ n 2 k\j2vri) barri- 
ers this reproduces Eq. ([S]) for delta barriers, y a , c — ct\.i- 
Analytical expressions for y a ^ c accurate to higher order 
can be easily obtained from the equation, which can then 
used to get the well depths in Eq. (|25]) . 

Likewise the amplitude of the particle to be inside the 
middle well is 



Ki sinh (kiD) 



(27) 
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which, for a narrow barrier n\D <C 1, also reduces to 
b w a7ro;i, the result of the delta barriers. 

Although the specific relations may vary with the size 
and shape of the potentials, the essential fact remains 
that the well depths need to be varied in sync with the 
barrier heights to maintain negligible population in the 
central well. It is interesting to note that by increasing 
the barrier widths and lowering its magnitude, keeping 
ai.2 constant, the probability of the particle to be in the 
middle well can be reduced further than for delta barriers, 
since the sinh function increases faster than linearly for 
larger k\D. 



IV. TRANSPORT PROCEDURE 

We apply the conditions derived above to demonstrate 
optimal transfer. We start with a particle in the ground 
state of well a which is isolated by a very large barrier 
from the neighboring well b (i.e., U\ — > oo, ol\ w 0). 

1. The barrier between b and c is lowered, i.e., «2 
reaches a finite (but small) value whereas o.\ w 0. 
At the same time the potential of well a is lowered 
to u a — —2/ct2, according to Eq. (|20)) . The state of 
the particle then coincides with the dark state hav- 
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FIG. 2: Numerical simulation of the time dependent proba- 
bility density for the particle confined to three wells. The pa- 
rameter l/ai,2 vary between 30 and 1550 as shown in Fig. [3] 
The plots correspond to times (a) t = 1400r, (b) t = 4, 900r, 
(c) t = 5,880r, (d) t = 7,000r, (e) t = 7, 700r, and (/) 
t = 14, OOOr with r = mL 2 /(ji 2 K) . 



ing Ad ~ —1, and Bd,Cd ~ 0. The two remaining 
states form superpositions of the states occupying 
wells b and c: B w C ~ l/\/2 (ground state), and 
Bw-Cwl/v^ (excited state). The energies of 
these two states are separated by ~ % ti^ailm from 
the energy of the middle state. So far the particle 
stays in well a. 

2. The barrier between a and b is gradually lowered to 
the same value as that between b and c, ot\ = 
with simultaneous lowering of the potential in well 
c, u c = u a = — 2/cv 2 - The state of the particle 
adiabatically passes to the superposition of occu- 
pying the wells a and c with the same probability. 
At the end of this step the particle is in the state 
with coefficients Ad ~ —Cd ~ 1/V% the remaining 
two states being separated from it by the energy 
of w spiU y^ax^lm. The probability of occupying 
the middle well b is in this stage ss 7r 2 af 2 /2- 

3. The barrier between b and c is gradually increased 
(0L2 — + 0) simultaneously bringing the potential en- 
ergy of well a to u a = 0. The middle state adia- 
batically goes to the state localized in well c, with 
coefficients Ad,Bd ~ 0, Cd ~ 1. The remaining 
two states now form superpositions with A«Brj 
1/V2 (ground state), and A w —B w -1/V2 (ex- 
cited state). The energies of these two states are 
separated by ~ h}k\a.\jra from the energy of the 
middle state. 

4. The barrier between a and b is increased to its orig- 
inal value, ai — > while bringing the potential en- 
ergy of well c back to u c = 0. Nothing happens to 
the particle which is now localized in well c. 

This procedure was implemented by a direct numerical 
solution of the time dependent Schrodinger equation by 
a split-step operator method, where the time propaga- 
tion is done alternately in incremental steps in position 
and momentum spaces, transforming between the two 
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FIG. 3: Time evolution of the parameters l/ai,2 characteriz- 
ing the magnitude of the delta potentials (see Eq. ([2|) with 
the time scale t = ml? /(n 2 Ti). The marks (a) through (/) 
correspond to the snapshots in Fig. [21 
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via Fast Fourier Transforms. Narrow Gaussian poten- 
tials were used instead of delta functions, based upon 
the considerations of the previous section. Snapshots of 
the evolution are shown in Fig. [5] and the correspond- 
ing variation of the barrier heights shown in Fig. [3] The 
simulation shows a smooth transition of the wavepackct 
from well a into well c, with essentially no population in 
the central-well at any time showing the effectiveness of 
the optimal conditions we derived. 



It is interesting to interpret the vanishing probability 
in the central well b in terms of the particle velocity. 
Classically, the current density is related to density and 
mean velocity by j — pv, from which the classical velocity 
in the middle of well b would be (assuming the simple case 
ai « c*2 and a'± w — a' 2 ) 



jb_ L<^2_ = 0^2 



(33) 



V. CURRENT IN THE CENTRAL WELL 

When the potential barriers are changing at a finite 
rate the wavefunction differs from the eigenstate of the 
instantaneous Hamiltonian. This change influences the 
resulting probability of occupying the middle well. The 
middle well occupancy and associated current can be es- 
timated using the continuity equation 



do 

Of 



.Hi. 

dx ' 



relating the probability density g = \ip\ 2 to the current 
density j. Integrating the probability density in wells a 
and c and assuming no accumulation of probability in 
well b it follows for the current j\ through well b to be 



3b 



d\c\ 2 d\a\ 2 2a.ia 2 (a 2 a\ ~ 
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(29) 



where the dot denotes time derivative. To allow for such 
a current in well b the sine-form wavefunction has to be 
complemented with an imaginary cosine part such that 



b\l — sin(fcox) 



ib*. 



— cos(fco£), (30) 



where 



ml? {a^di — aid 2 ) 
ir 2 h (a 2 + a 2 ) 3 / 2 ' 



(31) 



A dimensionless time derivative a' 12 = i~ai t 2 is defined 
in terms of the time scale r = mL 2 / (tt 2 ?i), which is of 
the order of the round-trip time T = 4mL 2 / (SttH) of the 
particle between the walls of an infinite well of length L: 



UiOIy — OL\tt 2 



a2)3/ 2 



(32) 



This value should satisfy the condition b% <C 1. In the 
case of u\ ~ a 2 and a\ « — a' 2 this leads to the require- 
ment a[ 2 <C a\ 2 - The dynamical contribution to the 
population of well b is less than the static one, I&2I 2 < \b\ 2 , 



provided that a' l2 < -Ka\ 2 . Even though for different 
time dependencies of ai^ 2 {t) we get different probabilities 
of occupying the middle well, we can use these results to 
estimate the trade-off between the highest probability pj, 
of being in well b and the time T necessary to complete 
the transport between a and c, namely pbT > nr. 



where vo = hko/m is the ground state velocity. Thus 
for a[ 2 Caj 2 (i.e., when the dynamical contribution to 
the middle-well probability is much less than the static 
one) the negligible probability of the particle to be in 
the central well can be understood as arising from the 
particle speeding through the middle well much faster 
than the velocity it has inside wells a and c. 



VI. CONCLUSIONS 

We presented analytical results for coherent tunneling 
via adiabatic passage in a system of three square wells 
separated by barriers that are high and narrow compared 
to the well dimensions. We have shown that in order 
to maintain negligible population in the central well, re- 
quired by the CTAP process, the depths of the extreme 
wells need to be varied as well as the barrier heights; this 
has the counterintuitive behavior that, during the time 
evolution, the depth of each exterior well needs to be cor- 
related with the height of the barrier non- adjacent to it. 
We determined those correlations along with expressions 
for the relevant stationary states and their energies sat- 
isfying them, which we apply in a numerical solution of 
the time dependent Schrodinger equation to demonstrate 
that they indeed lead to transfer between the extreme 
wells with virtually no occupation of the central well. 

Our general conclusion, that simultaneous and syn- 
chronized variation of both well depths and barrier 
heights is essential for optimal transfer, applies quali- 
tatively to more general potentials, although the exact 
nature of the relations may be different. For finite rate 
of time variation we found a relation for the tradeoff be- 
tween central well occupancy and the transfer rate. We 
also provide a novel interpretation of the transfer mech- 
anism based upon a current through the central well. 

The results here can find applications in coherent 
transport of electrons between quantum dots or atoms 
in micromagnetic traps, where our expressions can serve 
as a guide for choosing parameters for the CTAP process. 



Acknowledgments 

T.O. is supported by the Czech Institutional research 
plan MSM6198959213. K. D. acknowledges support of 
the Research Corporation in the initial stages of the work. 



7 



[1] J. K. Rowling, Harry Potter and the Order of the 

Phoenix, chap. 4 (Arthur A. Levine Books, 2003). 
[2] K. Bergmann, H. Theuer, and B. Shore, Rev. Mod. Phys. 

70, 1003 (1998). 
[3] K. Eckcrt, M. Lewenstein, R. Corbalan, G. Birkl, W. Ert- 

mer, and J. Mompart, Phys. Rev. A 70, 023606 (2004). 
[4] A. D. Greentree, J. H. Cole, A. R. Hamilton, L. C. L. 

Hollenberg, Phys. Rev. B 70, 235317 (2004). 
[5] J. H. Cole, A. D. Greentree, L. C. L. Hollenberg, and S. 

Das Sarma, Phys. Rev. B 77, 235418 (2008). 
[6] E. M. Graefe, H. J. Korsch, and D. Witthaut, Phys. Rev. 

A 73, 013617 (2006). 
[7] M. Rab, J. H. Cole, N. G. Parker, A. D. Greentree, L. 

C. L. Hollenberg, and A. M. Martin, Phys. Rev. A 77, 



061602R (2008). 
[8] J. Siewert, T. Brandes, and G. Falci, Opt. Commun. Opt. 

Commun. 264, 435 (2006). 
[9] L. C. L. Hollenberg, A. D. Greentree, A. G. Fowler, and 

C. J. Wellard, Phys. Rev. B 74, 045311 (2006). 
[10] S. Longhi, G. Delia Valle, M. Ornigotti, and P. Laporta, 

Phys. Rev. B 76, 201101 (2007). 
[11] W. Hansel, P. Hommelhoff, T. W. Hansch, and J. Reichel, 

Nature 413, 498 (2001). 
[12] D. Ferry and S. M. Goodnick Transport in Nanostruc- 

tures, Cambridge University Press (1997). 
[13] V.O. Nesterenko, A.N. Novikov, F.F. de Souza Cruz, and 

E.L. Lapolli, cond-mat/0809.5012 (2008). 



